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1 Introduction 



About 20 years ago the work on quantization of integrable models of Quantum Field Theory started with 
the idea of quantizing the classical soliton solutions [Q, 0|. Important results were achieved in this way, in 
particular for the sine-Gordon (SG) theory, the semi-classical spectrum of excitations (which happens to 
be exact quantum-mechanically) and the semi-classical approximation for the soliton S-matrix were found. 
The semi-classical S-matrix allowed to guess the exact S-matrix in the reflectionless case 0], and this was 
used further as a fundamental input in the bootstrap construction of the S-matrix for arbitrary coupling . 
Later, however the idea of direct quantization of solitons was abandoned in favor of other approaches such 
as Bethe Ansatz and its algebraic formulation in Quantum Inverse Scattering Method (QISM). 

Whatever the original motivations and methods were, it is fair to say that the most significant results 
in the theory were obtained by bootstrap methods. Exact S-matrices ||, U and exact form factors || [| 
were found by this method. Since the bootstrap calculation of the S-matrix is based on semi-classical results 
from quantization of solitons, and since the S-matrix defines the form factors through the set of form factor 
axioms ||, one may imagine that there is a direct path from quantization of solitons to form factors. Our 
attraction by this way of thinking does not mean, of course, that we reject the achievements of QISM which 
revealed the mathematical structure hidden behind the integrable models and led to discovery of quantum 
groups. 

In recent years two important pieces of information were added to the theory of sine-Gordon. The first 
of them is the relation with the perturbations of the minimal models of conformal field theory (CFT)j7[ B. 
The second is the discovery of the algebra of non-local charges, which is isomorphic to the (/-deformation of 
the universal enveloping algebra of the loop algebra s?2 H. Actually these two features are closely related: 



the restricted sine-Gordon (RSG) model coincides with the ^i^-perturbation of minimal CFT [1C|, 11 , but 
the restriction is intimately connected with the existence of the non-local symmetry. 

In the present paper we shall show that the results obtained by the bootstrap methods for sine-Gordon can 
be understood directly by quantizing solitons. Namely, we shall interpret the form factors of the restricted 
sine-Gordon model as matrix elements in a quantum mechanical n-soliton system Q This will allow us to 
underline the connection between profound structures in the classical and quantum theories: r- functions and 
separation of variables on the one hand and the space of the local fields on the other hand. Although we 
present the general structure for generic values of the coupling constant, we will reconstruct the sine-Gordon 
form factors only at the reflectionless points. We hope to return to the general case in another publication. 

For each n-soliton solution we shall introduce pairs of conjugated variables A, and Pi (i = l---,n), 
which in the quantum case satisfy Weyl commutation relations. Every local operator O can be considered as 
acting in this A-represcntation, and therefore can be identified with a certain operator 0(A, P). The typical 
formula for the matrix element of O between two n-soliton states can be presented as 

(B'\0\B) = J ^(A,B'yO(A,P)^(A,B) d^{A), (1) 

where ^(A, B) is the wave-function of the state of n solitons with momenta B%,- ■ ■ , B n . The measure d(i(A) 
will include a specific weight admitting a natural interpretation in the n-soliton symplectic geometry. We 
shall give explicit expressions for 0(A, P) corresponding to the Virasoro primary fields. 

In formula ([!]), the variables A will be complex. Choosing the integration domain is a non- trivial issue. 
It corresponds to choosing a real subvariety, which specifies the configuration space of the theory. This con- 
figuration space has a natural interpretation in the restricted sine-Gordon model where it can be understood 
as an analytical continuation of the sine-Gordon or KdV configuration space. The hermitian conjugation f is 
not the naive complex conjugation inherited from the sine-Gordon dynamics, but a more subtle one adapted 
to this choice of real subvariety. 

The local integrals of motion will be rewritten in the ^-representation. They are difference operators 
which can be expressed in terms of "quantum r-functions" . We will describe how to separate the variables 
in the associated Schrodinger equations. The fact that these Schrodinger equations are difference equations 
provides just enough room for the existence of the non-local charges commuting with the Hamiltonians. 
Recall that these charges do not allow direct classical limit. The importance of coexistence of two commuting 
Weyl subalgebras was pointed out in 

1 A side motivation for this study was to learn how to directly quantize solitons with potential applications to theories where 
soliton-like solutions (monopoles) are known. 
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2 The classical sine-Gordon theory. 



2.1 Sine-Gordon solitons. 

In this section we introduce a few useful notations for the sine-Gordon (SG) equation and its solutions. Let 
x± — x ± t be the light cone coordinates and d± — \{d x ±dt)- The sine-Gordon equation is: 

d+d-tp = 2sin(2<^) (2) 

It is convenient to introduce two r -functions t±, in terms of which the SG equations can be rewritten in 
Hirota form. The sine-Gordon field if is related to the r -functions by 

exp(i(f) = — 

r+ 

Let us describe the t -functions of the n-soliton solutions of the SG equation. Consider the function 

r{X 1 ,---,X n \B 1 ,---,B n ) = det(l + V) (3) 
where V is a n x n matrix with elements: 

V tj = 2 

Bi + Bj 

The n-soliton r-functions t±(x_, are written in terms of r as follows 

t±(x-,x + ) = T ± (X(X-,X+)\B) 

where 

t ± (X\B) = t(±X\B) 

The x±-dependence of X is quite simple: 

Xi{x+,x-) = X, exp(2(5^_ + Br 1 ^)) (4) 

The quantities Xi and Bi are the parameters of the solitons: f3i = log(Bi) are the rapidities and Xi are 
related to the positions. For the sine-Gordon equation, they satisfy specific reality conditions. For solitons 
or antisolitons, the rapidity B is real and X is purely imaginary, i.e. X = iee 1 with e = +1 for a soliton and 
e = —1 for an antisoliton. Wc shall not consider "breathers" in this paper but for completeness it should 
be mentioned that they correspond to pairs of complex conjugated rapidities (B, B) and positions (X, —X). 
Notice that these conditions are preserved by the dynamics. 

The sine-Gordon equation is a Hamiltonian system. The symplectic form is the canonical one: 



Visa — / dx Stt(x) A Sip(x) 

J — oo 

with ir(x) the momentum conjugated to the field f{x). Above, S denotes the variation on the phase space. 
The space of n-solitons can be viewed as a 2n-dimensional manifold embedded into the infinite-dimensional 
phase space. The restriction of the symplectic form on this finite-dimensional submanifold gives the n-soliton 
symplectic form. In the coordinates Xi and Bi it reads, cf eg.(l|, [l4| : 



\ ^ 



dX l dB, x - ABiBj dB t dBj 



A 



E^JDiJDj UDi UJDj 

2—1 2<J L J J 

The commuting conserved quantities are precisely the B's. In the following we shall always assume that the 
B have been ordered: < B\ < ■ ■ ■ < B n . A complete set of commuting Hamiltonians Hk can be chosen as 
the set of elementary symmetric functions 

H k =a k (B) (6) 

We recall that the symmetric functions crfc(B) are defined by: Yij( z + ^j) = J2k zn ~ k(J k{B). The local 
integrals of motion which are given on n-soliton solutions by 

n 

ii=s^ k+1) {B)^ B 7 {2k+x) 

3=1 
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can be expressed in terms of Hi. In particular for the light cone components of the energy- momentum we 
have J_ = If = Hi, and 1+ = If = H~ 1 H n -i. 

The variable Yj canonically conjugated to Bj is defined by Yj = Xj Tik^j ( a ' +bI ) ■ ^he symplectic form 

is then written as w = V) ■ A The equations of motion are very simple in the variables {Y, B\, which 

J I j Dj 

furthermore have the nice property of being separated. Other set of variables have also been introduced, in 



particular those which lead to Ruijsenaars models [115, 14 . However, for the purpose of comparison with the 



existing exact form factor formulae we need to introduce in the next section still another set of variables. 
A simple explicit expression for the r -functions can be obtained by expanding the determinant: 



t(x\b) = i + E II -\l x > (7) 

p=l rc{i, •■•■«} i<jel iel 

\I[=P 

= L + E^ + E ifUl.X.X, • ••• 

3 i<3 

with (3ij(B) = g' + g J . In appendix B, we gather a few useful formulae concerning these r-functions. In 
particular a very useful formula is the recursion relation satisfied by the n-solitons r-functions: 

t^(X\B) = T ^- 1 \X\B)+ T ^{0l l (B)X k \B) X n . (8) 
2.2 The analytical variables. 

We now give a parametrization of the n-soliton phase space in terms of new variables {^4, B}. The variables 
Bj can be considered as poles of the Jost function for the auxiliary linear problem, and the variables Aj 
are zeroes of the Jost function. The importance of these variables is better understood in the more general 
situation of quasi-periodic finite-zone solutions from which the soliton solutions are obtained by a limiting 
procedure. In the finite-zone case the analogues of Bj describe the moduli of the hyper-elliptic spectral 
curve, while Aj give the divisor of zeroes of the Baker- Akhiezer function. The general rule that the zeroes 
of the Baker-Akhiezer function give the correct set of variables for quantization was called the "magic 
prescription" in ]l6fl . We discuss the relation to the finite-zone solutions in Appendix A. Because of the nice 
algebro-geometrical meaning of these new variables we call them the analytical variables. 

The set of analytical variables {A, B} is related to the variables {X, B} introduced above by 

Bj - B k \ _ tt f Bj - At 

s B,+B k J J-J- 
k=i 3 yJ k/ k=i 



This relation can be considered as a system of equations for the symmetric functions <Jk{A) as functions of 
the {X, B} variables. The solution to this system is given in eq.(jsol) in Appendix B. 
In these analytical variables the symplectic form becomes 

We shall need the Liouville measure u> n = det(oj)dAi A • • ■ A dA n A dB\ A ■ • • A dB n with 

It is useful to write the non vanishing Poisson brackets 

r a o i _ U.kM B l Z Q^i^i ~ B ^ ( A2 r 2n 

Remark that the products in the right hand side can be written in terms of cross-ratios. This is the first 
manifestation of the conformal properties of these variables that will reappear in the following. 
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One can express the r-functions in terms of the variables Aj and Bj . The result is the following surpris- 
ingly compact formula 



The symplectic form as well as the r-functions enjoy an intriguing A <-> B duality. The proof of these 
formulae is given in the Appendix B. They lead to a formula expressing the sine-Gordon field in the {A, B} 
variables: 



r + XX y B 

Let us now introduce the variables Pj conjugated to Aj 



In terms of Aj and Pj the symplectic form takes the canonical form: 

„ dP,, dAj 
- = 2 ^lf A -Ai 

3=1 3 3 

We can express the hamiltonians Hk — Uk{B) in terms of the variables {A, P} using eqs.(|l4|) as a linear 
system for the Ok (B) . Surprisingly, the solution of these equations can be written in terms of the r- function 
as follows: 

The functions Tk are defined by 



r(Z\A) 3 ^ > 3 ll\ Aj -A k 



T k (Z 1 ,---,Z n \A 1 ,---,A n )= (16) 
J2 A n A l2 ■ ■ ■ Ai k T[Zi, • • • , —Zi t , • • • , — Z ik , • • • , Z n \Ai, • • • , A n ) 



il<i2<---<ik 



In particular t(Z\A) = t (Z\A) and T-(Z\A) = (U A j) r n (Z\A). We delay the proof of eq. @ as it turns 
out to be a limiting case of a more general quantum formula which we shall prove in Appendix E. It is 
convenient to introduce the generating function of the Tk- Let T(u) — J2k=o ukl ~k- It can again be expressed 
in terms of r-functions as: 

n 1/1 

T(u) = t(Z(u)\A) + uAj) W1 th Zj(u) = —^-L Z 3 (17) 



This follows from the quantum relation eq.(53) proved in Appendix E. 

Notice the unusual feature of our approach: the local integrals of motion ( |l5| ) are given in terms of 
r-functions. 

The equations of motion written in the variables Aj are as follows: 



d_A t = {I_,AA = H(A1 - Bj) J] * 



o + A t = {/+, = n itji n (is) 

j 3 m 1 3 
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These equations provide a particular case of the general equations of motion for the divisor of the zeroes of 
the Baker- Akhiezer function (see Appendix A). This kind of equations is commonly used in the theory of 
quasi-periodic solutions of integrable equations; they go back to Neumann and Kowalevskaya. One can show 
directly that ( |l8| ) together with ( |l3| ) imply the sine-Gordon equation ]l7], |l8| . 

To finish the discussion of the variables A4 we have to explain their trajectories. For one soliton the 
condition \e ltfi \ — 1 shows that the A-trajectory lies on the circle of radius B. Under the classical SG 
dynamics of n-soliton solutions every variable A/~ runs around a curve going in the lower half-plane from 
— Bk to Bk- When B\ <C B% <C • • • <C B n the trajectories are semi-circles, for finite Bk they are getting 
deformed, but not too much. For antisolitons these trajectories are replaced by their complex conjugate. 

2.3 Reduced action and the relation to the KdV equation 

Our goal is to quantize of the (R)SG theory in the soliton variables. The first step could be a semi-classical 
quantization. To perform it for n-soliton solutions, we need to compute the reduced action J q pdq since we 
are restricting the system to the level of Hamiltonians, as in the Maupertuis principle. The symplectic form 
( pi] ) can be rewritten as 

We discuss the relation of this symplectic form to the general theory of analytical Poisson structures |l9], 
in the Appendix A. The 1-form a is defined up to an exact form dF(A,B). In this paper, we shall only 
consider the possibility of a function F(A, B) independent of B (in order that a expends only on dAj, the 
coordinates) and of the special form F(Ai, ■ • • , A n ) — J^k Fk(Ak) (to preserve the separability property of 
the Aj). Hence the most general form of a is 

k=l \ j I k 

The functions Fk have to be fixed by additional considerations. 

Now we face a serious challenge: with these choices, the form a can not arise in the full SG theory for 
the simple reason that the reduced action S(Ai, • • • , A n ) constructed from a 

8<M> ■ 40 = t f k ^ (n f^) ■ t + Fk{Ak) 

can not be made real along the sine-Gordon soliton A, -trajectories by any choice of F^. This circumstance 
looks very discouraging and the variables Aj seem to be useless. But, as we shall later see, these are exactly 
the variables in which the comparison with the quantum form factor formulae is straightforward. We shall 
argue that the choices made for a correspond in fact to the restricted sine-Gordon (RSG) theory. For this, 
we first need to describe the connection with the KdV equation. 

Let us discuss briefly the relation between SG and KdV solitons. The KdV equation allows soliton 
solutions in the form 

u{x-) = d 2 _ logr(a;_,0) 

where r is exactly the same as before. In other words, the KdV soliton r-functions are identical to the SG 
soliton r-functions but with all the chiral coordinates x+ set to zero. The difference between the two cases 
lies in the reality conditions: the variables Xj which were imaginary for sine-Gordon become real in the KdV 
case. Before going into the details of these reality conditions, let us discuss the relation between the two 
cases in the fully complexified situation. 

In the variables Aj the equations for KdV solitons are exactly the same as for one chirality of SG (see 
Appendix A): 

e_4 = n(4 a - b j)U (20) 

3 j^i 1 j 
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The field u reads in these variables as 



K 3 3 ) 

It should be stressed that all the KdV fields are expressed in terms of even powers of Aj. The relation 
between the SG field ip and the KdV field u is given by Miura transformation: 

u(sc_) = -{d-(p{x-,0)f - id 2 _y{x-,<S) (22) 

It is a nice exercise to check directly that (^Tj) follows from e lip = JJ . -g 1 using the equations (p(f). 

As it has been said we are rather interested in the RSG than in the SG theory. The RSG model describe 
the <f>i.3-perturbation of the minimal model of conformal field theory (CFT). The reason why KdV equation 
appears in the context of two-dimensional CFT is well known: the second Poisson structure of KdV is a 
classical limit of the operator product expansion for the light-cone component of the energy-momentum 
tensor in CFT. (For the minimal models of CFT the classical limit is understood as the limit c — + — oo.) The 
second Poisson structure for KdV is: 

{u{x-),u(x'_)} = 6'(x--x'_)(u(x-)+u(x'_))+5"'(x--x'_) (23) 

The KdV field u is identified with the classical limit of T As explained in Appendix A, it turns out that 

the second KdV Poisson structure restricted to the soliton manifold in the {A, B} variables is identical to 
the symplectic structure ( |Io| ) that we have derived from the sine-Gordon theory. It is remarkable that the 
conformal KdV Poisson structure appears in the (massive) sine-Gordon model when restricted to the soliton 
solutions. Again, this is only true on the soliton submanifold. 

Thus we are considering the light-cone hamiltonian picture, and the lines x + = const are the space 
directions. The coordinate x+ has to be considered as time. Globally we can not introduce the x+ dynamics 
in KdV theory, but this can perfectly be done on the n-soliton solutions: the corresponding Hamiltonian is 
I + = ")2 Bj l and the equations of motion have the familiar form 

^=n^^n^ (24) 

3 jYi • 3 

For real solutions the x+-dynamics of Aj is organized as follows (see Appendix A). Each Aj (i < n — 1) 
moves inside either of two segments (Bi,Bi + i) and (— -Bj+i, — -Bj), the points Bj and — Bj being identi- 
fied (recall that the observables depend on A\ only). The point A n moves in two segments (B n ,oo) and 
(— oo, —B n ) the points ±oo being identified, notice that oo is a regular point for the equations (pi[). 

Let us now normalize the 1-form ([19]) in order that it is real over x+-trajectories. To do that we have first 
to fix the branch of logarithms. For further convenience we put n cuts over semi-circles in the upper-half 
plane connecting the points Bj and —Bj and require that the logarithm 

is real when — B\ < A < B\. With these conventions one easily figures out that the following choice of Fk 
corresponds to a 1-form a, real over ^-trajectories: 

. dAu 

kni | • -± (25) 

The reduced action is real as well because the trajectories are real. 

Returning to the relation with the sine-Gordon model, it is obvious from ( p2] ) that real SG solutions do 
not correspond to real KdV solutions. One important consequence of the comparison with the form factor 
formula is that we are actually considering an analytical continuation of the sine-Gordon and KdV dynamics. 
More precisely, the semi-classical quantization exactly brings out the 1-form a as defined above in eq.(Eq). 
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But the trajectories for A^ lies on semi-circle of radius Bk < \Ak\ < Bk+i- In other words, the 1-form 
needed for the quantization of RSG model is obtained by analytical continuation of the 1-form (p5|). In this 
analytical continuation we identify the variable Ai which runs from Bi to — Bi in SG case with the variable 
running inside {Bi, -Bi+i) U {—Bi, — 2?j +1 ) in KdV case (we put B n+ i = oo). Selecting the trajectories in the 
complex A-plane is just choosing the real phase space of the theory. 

Following our logic the form ( p5| ) must be understood as restriction to n-soliton submanifold of a globally 
defined 1-form related to the Poisson structure (p3|). We are not able to define such a global 1-form, but 
we conjecture that it exists. We can only refer to certain self-consistency checks (such as reality of ( ]25| ) ) to 
support this conjecture. 



3 The SG form factors in absence of reflection of solitons. 

3.1 Sine-Gordon versus Restricted Sine-Gordon. 

The sine-Gordon equation follows from the action: 



S=- £ dfx, C = {d.ifY + m^(cos(2<p) - 1) 
1 J 

where 7 is the coupling constant, < 7 < 7T. The free fermion point is at 7 = ^. In the quantum theory, 
the relevant coupling constant is: 

= 7T7 
7T — 7 

We shall always use the constant £, which plays the role of the Planck constant. Only the mass is renormalized 
but not the coupling constant 7 

The SG theory is invariant under the quantum affine loop algebra t/^s^) with q = cxp(i^|-). This 
symmetry algebra is generated by the topological charge and four non-local charges of spin s — i^. (By 
convention, q here is the inverse square of that of ref.§.) 

The canonical stress tensor of the SG theory is T~„ = {d ll ip){d^(p) — ^-C. The SG theory contains 
two subalgebras of local operators which, as operator algebras are generated by cxp(«<^) and exp(— itp) 
respectively. Let us concentrate on one of them, say the one generated by exp{iip). It is known that this 
subalgebra can be considered independently of the rest of the operators as the operator algebra of the theory 
with the modified energy-momentum tensor: 



G 



Tp a = T$T + in^j ^—^e^e^d^ip (26) 

This modification changes the trace of the stress tensor which is now: 

T ™^ = m 2 exp(2^) 

This modification corresponds to the restricted Sine-Gordon theory (RSG). For rational - it describes the 
$i 3-perturbations of the minimal models of CFT, but it can be considered for generic values of £ as well. 
The central charge and the dimension of the operator $13 are given by 

1 671-2 a £- 7r 



Modifying the stress tensor as in eq. ( p6| ) modifies the Lorentz boost and hence the spin of the non-local 
charges. Under this modification two of these charges become spinless. Together with the topological charge 
they then form a representation of the finite quantum algebra U-^sh)- The restriction of the Sine-Gordon 
model consists in gauging away this symmetry subalgebra. The physical states of the RSG model are 
annihilated by these spin-less non-local charges. The physical operators of the RSG model are those which 
commute with these charges. In particular, e lv which commute with them is a physical operator, but e~ llf 
is not since it does not commute with these charges. 

The asymptotic states of the RSG theory are constructed as follows. Consider the states in the SG theory 
containing n solitons and n anti-solitons: 

|/3l,/?2, • • • ,/32n)e 1 ,e 2 ---,e 2 „ 
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where 0j are the rapidities of particles, and tj is + or — for soliton or anti-soliton respectively. For the RSG 
theory one introduces the states 

\Pl,02, ' ' ■ ,02n))e 1 ,e 2 -,e 2n = eX P(^| X] ' ' ' , #2ra) ei ,e 2 -,e 2 „ (27) 



The extra factor in \0\,02, • ■ ■ ,02n}) is an echo of the modification (26) of the stress tensor. In the new basis 
the S-matrix becomes manifestly invariant under t/^s^)- The asymptotic states of the RSG model are then 
the L^(s/2)-scalar in the Hilbert space span by the states \0i, 02, • ■ ■ , fan))- 

In the restricted theories one generally can not introduce a positively defined Hermitian structure. Obvi- 
ously the SG hermitian conjugation * maps the RSG model into the symmetric restricted model constructed 
from exp(— tip). On the other hand the SG charge conjugation (ip — > — tp) is also broken in the RSG model. 
However, one can introduce an anti-linear involution for the RSG model as the combined CT reflection. For 
any local operator the combined transformation (which we denote by f) corresponds to hermitian conjuga- 
tion and reflection tp — > —tp. This operation f does not give a positively defined scalar product for the SG 
theory, but it does not lead to contradiction if one stays within the RSG theory. In particular, since the 
modified stress tensor T™° d is hermitic in the RSG theory, one has: 

This is a simple but fundamental remark. 
3.2 Form factor formulae. 

In what follows we shall consider the case £ = ~ for v = 1,2, when the reflection of solitons and 
anti-solitons is absent. In these cases, the two parameters q and q are 

q = e^ , q = l 

The S-matrix is diagonal and given by 

8(3) TT sinh ^ + ^ } - 'ff (^-^ 

We shall use the following notations: 

2tt 

B = exp(/3), b = exp(— 0) = exp(2^/3) 
Consider any local operator 0(x) for RSG. Its form factors are defined by 

fo(0U02, • • • , /3 2 nk,e 2 ..,e 2 „ = «0|O(0)|/?i, ft, ■ • • , 02n))er,e 2 -,e 2n (29) 

The form factors are given by the formulae 

xfo(0i,0 2 ,---,02n)-.- + ...+ (30) 

The function ((0) is regular in the physical strip < Im < ir, and satisfies £(— 0) = C(0—2in) = S(0)<^(0). 
It can be found in |j) . We shall not consider this prefactor in eq.(30) since it is the same for all operators: 
it is related to normalization of the wave function which we hope to discuss sometime. The most interesting 
part of the form factor is given by 



fo(01,02, ■ ■ ■ ,02n) + •■■ + 

-if p n 2n n 

—— j dAf- j dA^W^iUB^WiA* - A*) L (A 1 ,---,A n \B 1 ,---,B 2n )]lar i |.TI ) 

r, »=lj = l i<j i=l 
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Some comments are needed for this formula. The function ip(A, B) is given by 



^(A,B) = H(B-Aq-1) (32) 

3=1 



This function satisfies the difference equation: 

As usual we define a = A 2l/ . The contour C is drawn around the point A = 0. Notice that the right hand 

side of the formula (|3l|) does not depend on the partition of the particles into solitons and anti-solitons, this 

is a peculiarity of reflectionless case. 

Different local operators are defined by different functions Lo(Ai, • • ■ , A n \B\, • • • , B 2n ). These functions 

are symmetric polynomials of A\, • • • , A n . For the primary operators <fr 2 k = exp(2ki(p) and their Virasoro 

descendants, Lq are symmetric Laurent polynomials of B\, ■ ■ ■ , B 2n . For the primary operators <&2fc+i = 

i 

exp((2fc + they are symmetric Laurent polynomials of B\, • • • , B2 n multiplied by Y[ B? . Our definition 

of the fields <I> m is related to the notations coming from CFT as follows: <I> m corresponds to $i, m +i. The 
explicit form of the polynomials Lq for the primary operators is as follows 

n 2n 

L* m (A ir --,A n \B 1: ---, B 2n ) = H A? H B7* 

We explain in Appendix C that this definition agrees with the formulae for the form factors of the operators 
$i and $2 given in We shall return to this definition in the following sections. 

Usually the formulae for the form factors are written in a slightly different way || First, the 

integrals are (n — l)-fold, second, instead of the polynomials IK^f — A?)Lo{A\, • • • , A n \Bi, • • • , B 2n ) under 
the integral, one usually has polynomials of the type Y\{Ai — Aj)Lo{A\, • • • , A„_i|i?i, • • • , £?2„) with the 
limitation degAi{Lo) < n - The conventional form of the integrands is important for generalization to the 
case of arbitrary coupling constant because, in this case, the integrals become more complicated and the 
above limitation is needed for convergence. In Appendix C we explain briefly how to rewrite formula ( |3l| ) 
in the conventional way. 

There are two important facts that we learn from the calculations of Appendix C. 

1. The substitution into ( |3l| ) of the polynomial Li(A\, ■ ■ ■ , A n ) = 1 which corresponds to the unit 
operator, gives zero because there is no simple pole in the contour integral, in agreement with the fact that 
the matrix element of the unit operator between the vacuum and an excited state vanishes. Also if we try 
to substitute into (^l|) the functions 

L*_ m (At, ■ ■ ■ , A n \B u ■ ■ ■ , B 2n ) = J] Aj m Y[ BJ 

3 3 

which must correspond to the operators <5>- m = exp(—miip) the integral vanishes. This means exactly that 
our formulae suit rather the RSG than the SG model. This does not mean that the formulae for the form 
factors of <&- m do not exist: they are obtained by SG charge conjugation, we want to say only that these 
formulae can not be obtained by putting L$_ m into the integral formula. 

2. Consider any polynomial M(Ai; A 2 , ■ ■ ■ , A n ) anti-symmetric with respect to A 2 , - ■ ■ , A n . The value of 
the integral ( |3l| ) does not change if we add to LJ(^i — A^)Lq(Ai, • • • , A n ) an "exact form" 

l^iA^A!, ■■■..{,.■■■. An) Y[{Bj + A k ) - qM{qA k -A u ■■■,A k ,---,A n ) ]]_(B 3 - A k )^j (34) 

where A k means that Ah is omitted. 



4 The semi-classical analysis. 



10 



4.1 The semi-classical analytical quantization. 

In this section we show that the exact formulae for the form factors in the semi-classical limit can be obtained 
from the semi-classical quantization of solitons in the analytical variables. We shall proceed to the exact 
quantization of solitons in the next section. 

We shall consider the matrix elements of the operator O calculated between two n-soliton states, instead 
of the ones calculated between the vacuum and the state with n solitons and n anti-solitons. This is because 
they have a more direct semi-classical interpretation. Crossing symmetry relates the connected parts of these 
two kinds of matrix elements. The essential piece of the n solitons to n solitons form factor is given by 

fo(P[,---,(3' n \Pi,---,p n ) = ^-^ J dAf-j dA n J{'\{i>{A i ,-B' j )^A i ,B j )Jl{A 2 i -A])\{aT i 



c c 



i=lj=l i<j 



xL (A u - ■■,A n \-B' l ,---, -B' n , Bx, ■ ■ ■ , B n ) (35) 
where ip(A, B) is defined in eq.(|3^). We would like to present this formula as 

*(A B')^0(A, P)*(A, B)W{A)dA (36) 



where the various terms should be interpreted as follows. The function ^{A, B) is the n-soliton wave function 
in the Schrodinger ^-representation. The scalar product in the Hilbert space is written in terms of a non 
trivial weight jy(A) and a possibly non standard conjugation f. The integration domain T> has to be 
compared to the configuration space of the classical theory. The quantity 0(A, P) is the operator O realized 
in terms of A and the conjugated variables P. The detailed discussion of the operators A and P will be 
given in the next section. Let us first make a semi-classical step in this direction. 

Consider the semi-classical approximation of the function ip(A, B). Asymptotically we have (see 
Appendix D) 

i>(A,B) ~™ . exp ( - /log (^-4) A^dA I OiTi 



y/B 2 - A 2 °\B + A 

\ o 

where the integral is taken over a contour which does not cross the cut going along the semi-circle in the 
upper-half plane from B to —B. The logarithm is real when — B < A < B. An important property of the 
reflectionless situation is that the integrals over different contours of this kind give the same result when 
substituted into the exponent. Indeed we have, 

log ( 1-4") A-'dA = -27T 2 



7 



B + A 



where the contour 7 is drawn around the cut. Therefore, when v is integer there is no ambiguity in the 
exponential in eq.(p7|). This is similar to the topological ambiguity of the WZNW action. 

Let us now construct the semi-classical wave-function in the ^-representation for given values of B. The 
general rule is |2l]| : 

V(A) ~ {da(A)Y exp ^rS(A) (38) 
in 

where d[i(A) = uo 11 is the Liouville measure at the point A. Let us give some explanation. Integrable models 
provide examples of a situation which is usually described in textbooks only for the case of one degree of 
freedom (i.e. when there is only one integral of motion: the energy). In our case in order to construct the 
semi-classical wave-function in the A-representation, we have to consider the polarization corresponding to 
the A variables and the integrals of motion B. In these variables the Liouville measure can be written as 

d(i(A, B) = [dot uj(A, B)\ dA x A • • • A dA n A dB x A • • • A dB n 

So (d/j,(A))^ has to be understood as (detu>(A, B))i . In this measure the terms which depend only on B 
will be omitted since they correspond to the normalization of the wave-functions which we do not consider. 
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In eq.(|38|), S(A) is the semi-classical action. We choose as semi-classical action the analytical continuation 
of the KdV reduced action discussed in the previous section, which in the analytical variables A is S(A) = 
J2k s k(Ak) with 

Sk(A k ) = J log (j[ ■ ^ - ™klog(A k ) (39) 

where the integrals are again taken in the plane with cuts going along the semi-circle from Bj to — Bj for 
all j. The ambiguity in the definition of the contours is irrelevant for semi-classical quantization at the 
reflectionless points as it has been explained above. For the case of generic coupling constant we can not 
manage with the variables Ai, instead the variables a« = log(-Aj) will have to be considered. We shall return 
to this point in future publication. 

We now have all the necessary ingredients: the classical action and the Liouville measure. So, up to 
normalization depending only on the Bj , the semi classical wave function is 

%(A 1 ,--' ) A n \B 1 ,--- ) B n )~ (40) 

nn(A 2 - B])j cx P U (^t / log (n, • # - ^ma^ j 

The Planck constant is identified with £ = -. We shall divide ^ into the wave function ^ and a 
independent piece which will be put into the integration measure: 

*(A!, • • • , A n \B u •••,£„) = (W(Ai, • • • , A n ))H{A u ■ ■ -,A n \B u • • • , B n ) 

where 

n 

W{A x ,---,A n ) = \{{Aj - A^)\{a- k 

i<j k=l 
n n 

^(Ai,- ■ ■ , A n \Bi, • • • , B n ) = nn^>^) 

where we recall that a k — Ai". 

Notice that choosing the reduced action as we did, i.e. as the analytical continuation of the KdV action, 
was crucial for producing the factors ]J k k in the weight W^(A). Identifying ^{A, B)^ = ^(A,—B) and 
the integration domain T> with the contour integrals over C we see that the form factors in the semi classical 
limit admit a representation as in eq. (^6|) . 

It remains to find a more physical and geometrical interpretation of the integration domain and of the 
conjugation f in the semi-classical regime. Under the classical dynamics of n-soliton solutions, every variable 
Ak runs around a curve going in the lower half-plane anti-clockwise from — B/~ to B^. As we already know 
there are two difficulties: 

1. The action s(A k ) is not real along the classical trajectory. 

2. The trajectory itself is complex, and its exact geometry depends on B^ which must become the spectral 
data after the quantization. 

On the other hand the exact quantum formula show that the integration domains should be small circle 
around the origin. This suggests the following prescription for fixing these problems. 

1. Because of quantum-mechanically possible penetration through the barrier one has to consider the closed 
trajectories, i.e. to add to the classical trajectory of A^ the piece running in upper half-plane anti-clockwise 
from Bk to — B^. Classically this extra piece corresponds to an anti-soliton trajectory. 

2. In order to eliminate the geometrical dependence of the trajectory of Ak on Bj (leaving the topological 
one) we replace it by arbitrary closed curve along which Bk < \A k \ < B k+ x, There is one more reason for 
considering the configuration space as the one composed of closed cycles. If we regularize the theory putting 
it onto large but finite interval with periodic boundary conditions the trajectories of Aj will become the 
closed cycles of the same type as described above (see Appendix A). 

3. Last but not least, the reality condition is taken as follows: 



s k (A k ) = Sk {A k ) = s(A k y (41) 
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Let us explain that this condition is satisfied by (tS9j) . If we take 





over conjugated contours, then ( [i"l| ) would be true if the integrands were continuous over these contours. 
However, we defined the logarithms in such a way that they have cuts in the upper half plane, and to reach 



Ak we have to cross k cuts. The addition of — nik \og(Ak) in (39) is needed to compensate the corresponding 
discrepancy. 

The function vfft is constructed from s'. The choice of the conjugation f ensures that ty(A\B'y = 
#(A| - B') as it should be. 

Let us discuss the condition (f4l|). This condition provides the usual reality condition 



Sk(A k ) = s k (A k ) 

if Ak is real and such that B k < Ak < Bk+i- The latter segment is the trajectory of Ak which corresponds 
to real n-soliton solutions of KdV equation. So, we are dealing with an analytical continuation of KdV 
mechanics, as it has been explained in the Section 2.3. This fits perfectly with our assumption that we 
are actually doing not SG but RSG model, which has to be considered as a continuation of KdV. Also the 
fact that the involution ([ll]) connects the soliton and the anti-soliton part of the trajectory has a physical 
explanation: as we already said RSG is neither C nor T invariant, so one has to consider the combined CT 
invariance. The involution (|4l]) corresponds exactly to the combined CT reflection. 

It now is clear that the formula ( |35| ) allows a semi-classical interpretation by means of this semi-classical 
wave function. We stress once again that this semi-classical analysis forces us to introduce the analytical 
continuation (|d|). The insertion of a local operator corresponds to the insertion of a polynomial Lq which 
can also be understood from the expression of the classical local observables in n-soliton variables. We shall 
return to this point in the section on exact quantization. 

Let us emphasize one circumstance. In order to put together two semi-classical wave- functions ^(A, B) 
and ^(A, B'y into the matrix element, one has to be sure that the contours of Ak can be identified. Recall 
that for these two wave functions the contours are inside Bk < \A k \ < Bk+i and B' k < \A k \ < B' k+1 
respectively. This shows that the semi-classical quantization works nicely if the soliton states are not very 
far from each other in the following sense. The sets Bj and B'j are ordered Bk < B k+ i and B' k < B' k+l . The 
semi-classical formulae are applicable if in addition B' k < B k +i and Bk < B' k+1 . In the next section we shall 
discuss this point further. 



4.2 Stationary point calculation of the semi-classical integrals. 

In this section we shall show that the semi-classical integrals for the form factors are actually given by 
stationary point contributions. This result, which is interesting by itself, will also provide some justification 
of the compactification of the configuration space made before. 

It is convenient to return to the form factor taken between the vacuum and the state with n solitons and 
n anti-solitons. One can show that the semi-classical limit allows the necessary analytical continuation. Let 
us consider the exact quantum formula: 



fo(Pl,02, ■ ■ ■ ,/hn) + ■■■+ 

n 2n 



1 /* nib ^11 lb 

— JdA l ---J dA n l[l[^(A i) B j )l[(A 2 i -A 2 j )Lo(A 1 ,---,A n \B 1 ,--- > B 2n )l[a- 

n ^ r i=l j'=l i<j i=l 



(2m) 

Ci 



(42) 



In quantum formula the contours of integration are arbitrary. However, to perform the classical limit we need 
to take them as follows: if A k £ C k then B^k < \A k \ < i?2fc+i- The point is that the asymptotic of ip(A, B) is 
not a single- valued function of A, so one can not move the contours of integration after the asymptotic of the 
integrand is calculated. The above prescription for the contours agrees with the semi-classical considerations 
of the previous section. The justification of this choice of contours follows from further stationary phase 
calculations. The momentum Bi corresponds to either soliton or anti-soliton. The matrix element is related 
to that considered in the previous section by analytical continuation. The concluding remarks of the previous 
section show that it is exactly the choice of cycles described above which corresponds to the semi-classical 
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matrix element, and the momenta Bj are partitioned into pairs Bn-\,Bn such that in every pair we find 
one soliton and one anti-soliton. 

The asymptotic of the integrand is written down using the formulae from the previous section. Let us 
consider the integral with respect to Ak- It contains the divergent exponent: 

«p (h (/^n (frt) dAk - ^*^))) 

where the last term comes from 7J fc a^ k . One has to consider the stationary phase point at A = —iFk which 
solves the equation 




= 2nik 



It is easy to see that all Fk are real and positive, moreover B^ < Fk < -B2fe+i for k = 1, • • • n — 1, F n = oo. 
The calculation of the second derivatives shows that the integrals over Ak for k — 1, • • • n — 1 are given by 
stationary point contributions coming from — iFk if the contour of integration goes trough the point Fk being 
topologically equivalent to a circle lying inside the domain B^k < \Ak\ < £?2fe+i' Obviously the contours Ck 
from ( |42| ) can be drawn like this. The point oo for the n-th integral is not a true stationary point, but for 
the n-th integral the divergent exponent vanishes when A n — > oo, and the integral is sitting on the residue 
just like in the exact quantum calculation (Appendix C). 

The equations for the stationary points can be summarized into the Bethe Ansatz type equation 

Consider any polynomial M (A±; A2, ■ • • , A n ) anti-symmetric with respect to A2, ■ • • , A n . The main con- 
tribution into the asymptotic of the integral ( p^ ) does not change if we add to Y\{Af — Aj)Lo(Ai, ■ ■ ■ , A n ) 
the expression 

^2(-l) k M(A k ; At, • • • ,Ak ,■ ■ ■ , AOfll^ + Ak ) - 11^ - Ak )) ( 43 ) 

k 

because it vanishes on the stationary point if k < n — 1, and cancels residue if k = n. Comparing this 
formula with exact quantum formula ( |34| ) we observe the following nice circumstance. The semi-classics is 
of course not exact like it happens in geometrical quantization, however one can think of the exact quantum 
relation (Q) as of deformation of J43| ) which in turn is nothing but the consequence of the equation for the 
stationary points. The factorization over (^) and ( |34| ) leave the same number of independent polynomials 
of A k . 

The stationary phase calculation fits nicely with the classical picture. The point —iFk lies in the region 
acceptable by analytical continuation from the classical soliton trajectory. This fact gives an important 
justification of our semi-classical methods. 



5 Exact analytical quantization of solitons. 
5.1 Hilbert Space and Hermitian structure. 

In this section we do not specify £ to the values - until the last subsection in which we shall reconstruct the 
RSG form factors. Consider the variables otj = \og{Aj) and Wj = log(Pj). The symplectic form is canonical 
in terms of these variables: 

n 

lo = 2 d/Kj A dctj 
i=i 

We quantize them in a canonical way: oij act by multiplication and nj by differentiation, i.e. aj — * aj and 
TTj — > gf-- Here £ plays the role of Planck constant. The operators A and P are defined as: 



Aj = exp(ctj), P 3 = exp(i£-g— ) 



d_ 

h 
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They satisfy Weyl commutation relations with q = cxp(i£): 



Pj • Aj = qAj-Pj, with q = e li (44) 

Pk-Aj = Aj-Pk, j^k 

This definition is rather formal since we have not yet defined the Hilbert space of functions of aj. Staying 
on the same formal level one realizes the following important circumstance Jl2| . There is another pair of 
operators 

2tt . d . 

a, = exp(— aj), Pj = exp(27ri^— ) 

which represent the Weyl algebra but with the dual quantum parameter q = exp(^-^), associated to the 
non-local symmetry algebra U^sl^)'- 

~2 , - , 27r V 

Pj ■ aj — q aj -pj, with q — exp( - ) 

Pk-aj = aj-p k , j^k (45) 

The operators aj and pj commute with Aj and Pj . The existence of two dual algebras will be crucial in the 
following. 

All these operators act on wave functions $>(a). In the reflectionless case "J/ (a) is a single- valued function 
of A's so we shall write ^(A) in that case. The operators Aj act by multiplication whereas Pj act by shifting 
the argument of \& by i£: 

Pj*(ai, ■ • • ,aj,- ■ ■ ,a n ) = • • • ,aj ■ ■ ■ ,a n ) 

To complete the representation of the canonical commutation relations in the Hilbert space, we also need 
to introduce the scalar product. In order to take into account the reality condition (e lv y = e ltp ', specific to 
the RSG model, we define the scalar product in rather unusual way. Let ^>i(A) and ^(^l) be two wave 
functions, then: 

(*i|* 2 >= f day - f da n ]jA i ]l(A?-A?)y{(a)Q(ai,---an,Pi,---,Pn)Ma) (46) 

•' G ^ C i i<j 

where the hermitian conjugate wave function w|(a) is defined by: 

^\{a)=WJ^) (47) 

Notice that although this definition involves two complex conjugations, it is anti-linear as it should be. 
The operator Q{a\, ■ ■ ■ a n ,pi, •• • ,p n ) inserted into the integral is not specified here, except for a simple 
constraint arising from the condition | ^2) = (XE r 2 1 ^1) ■ The exact form of this operator has to be fixed 
from additional requirements. The role of this operator is similar to that of screening operators in CFT. 
This analogy is not just a formal coincidence as it will be explained later. Due to commutativity of A, P 
with a,p, the particular form of the polynomial Q is irrelevant for formal properties of the operators A, P. 
The contour C is complicated for generic £, but reduces to small circle around the origin in the reflectionless 
case. 

As usual, given an operator O, its hermitic conjugated operator is defined by ((C^i)|\E , 2) = 
I (C\I/2 )) • For the canonical operators Aj and Pj this gives: 

(48) 

^=<lPr ]| i "V I (49) 





It is an interesting check to verify that these relations are compatible with the Weyl commutation relations 
d44|). The relation for Aj is obvious from the definition. The formula for Pj can be deduced as follows: 

(*i|(P 3 -*a)) - 
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= / dai ■ ■ ■ I da n ]^[ A t - A ]) *i(a)Q(ai, • • -a ni pi, ■ ■ ■ ,p„)* 2 (- • • ,&j + if, ■ ■ •) 

^ C ^ C i i<j 

= qf dar ■ ■ ■ f da n ]jA t ]J(A^ - A]) J] ( * t j ~ £° ) • • ,a> + W, ■ • 



xQ(ai, ■ ■ ■ a n ,pi, ■ ■ ■ ,p„)* 2 (a) 

The second equality follows from the first by changing variables ay — > ay + if which is possible if \& is regular 
as function of ay between C and C + if. To obtain these formulae we crucially use the fact that q is of 
modulus one, ic. q~ x = q. (They will not be correct for real values of q.) 

The fact that Aj is hermitian simply expresses the fact that the field e l(p is a real field in the RSG theory 
since it has to be identified with the field $13 of the minimal conformal model. In other words, the quantum 
variables Aj are real variables leaving on the unit circle! 

5.2 Hamiltonians and quantum r-functions. 

We now introduce the quantum version of the hamiltonians defined in eq. (|l5|) . Since the classical 
hamiltonians are complicated functions of A and P, we have to specify the order of the operators in the 
quantum formula. We choose the following minimal deformation: 

T{P'\A) H k =%{P'\A) (50) 

with P' = (-l) n P and 

T k (P\A) = T{P 1 ,---,-P ll ,---,-P lkl --- 1 P n \A) A ll A l2 ---A lk (51) 

ii <%2 - <ik 

A il A i2 ---A ik T{P 1 ,---,-qP il ,---,-qP ik ,---,P n \A) 

il <%2---<ik 

We shall refer to the operator T{P\A) as the 'quantum r-function'. Its explicit expression is: 

n 

t(p\a) =1+^-1)^ e cn^) n p < ^ 

P=l fC{l,-,»} iei ig/ 

|/|=p 3<tl 

with 

qAj + Aj 



Ho 



Ai Aj 



Alternatively, T(P\A) can be recursively defined by eq. (p4|) in Appendix E. In the limit q — > 1 this formula 
coincides with the classical r-functions t(Z\A) as defined in eq.(0). For two particles: 

The operators T(P\A), but not the hamiltonians H n , are closely related to the Mac Donald difference oper- 
ators J23|. In particular the terms in T(P\A) of fixed homogeneity degree in P form a family of commuting 
difference operators. More precisely, the two generating functions T{uP\, • • • , uP n \A) and T{vP\, • • • , vP n \A) 
commute for any u and v. Note also that T n (P\A) = T(-P\A) ■ (11^) = Ai)T(-qP\A). It is conve- 
nient to introduce the generating function T(u) for the quantum r-functions T k : T(u) — XX=o uk ^k- As 
shown in Appendix E, we have the following expression for T(u): 

n 14 

f(u)=T(P(u)\A)-l[(l + uA j ) with />•(„) =P J .-^Li (53) 
3=1 3 

This can be used to relate the generating functions of the symmetric functions of the B and A operators as: 

n n 

+ uBj) ■ T(P\A) = T(P(u)\A) ■ JJ(1 + uA,) 
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It is important to realize that all the hamiltonians Hj~ commute with the operators Oj and pj , since they 
are functions of the A and P only. As we already pointed out, the operators Oj and pj are associated to the 
quantum afhne symmetry U-^sl-z). Thus, the fact that the Schrodinger equation is a difference equation just 
leave enough room for this non-local symmetry. This symmetry does not have any straightforward classical 
meaning, it corresponds to choice of topologically different components of the classical configuration space 
(different cycles). The way of encoding this topological information into quantum formulae through the 
algebra which commute with all local observables is, in our opinion, a very interesting feature of quantization 
of solitons. 

The conditions for the H k to be hermitian are bilinear identities on the quantum r-functions : 

T{P\A) T k (P\Ay = %{P\A) T{P\A)^ (54) 

Furthermore, the quantum r-functions behave nicely under hermitian conjugation. More precisely, the 
definition of pj implies that the quantum r-functions are not hermitian for our scalar product but that their 
hermitian conjugates are again quantum r-functions: 

T(A 1 P 1 , ■ • • , A„P„|A)t = TfaAiPi, • ■ ■ , qKP n \A) 

for any real parameters Ai, • • • , A n . Here, once again we use the fact that q is of modulus one. In particular, 
T(P\Ay — T(qP\A). This allows us to rewrite eq.(|54|) directly in terms of the quantum r-functions. For 
example, for k = n we have T n (P | A)t = Q1, Ai)^T(-P\A)^ = fl^ Ai)T(-qP\A). The hermiticity condition 
of H n is then equivalent to: 

T(P\A) T(-P\A)=T(-P\A) T{P\A) 

which is a consequence of the identification of the quantum tau function as the generating function of the 
Mac Donald difference operators. Besides the case k = n, we also checked the hermiticity conditions ( |5^ ) 
for all two-particle hamiltonians. 

The conditions that the hamiltonians Hk commute can also be rewritten as bilinear identities for the 
quantum r-functions: 

%{P\A) 7I(P|A)t = Tt(P\A) %{P\A)^ (55) 

To rewrite these conditions in a simple form, we used the hermiticity of the hamiltonians. The hermiticity 
condition ( |54| ) corresponds to eq.(|55|) with I — 0. Eq.(|55|) can be rewritten as a bilinear identity on the 
quantum generating function T(u), cf Appendix E, eq.(J8q). We are missing a complete algebraic proof of 
these conditions for arbitrary n. Although we don't have any doubt that it is true because the hamiltonians 
admit simultaneous eigenfunctions. 



5.3 Quantum separation of variables and soliton wave functions. 

One of the magic aspect of these quantum hamiltonians is that the corresponding Schrodinger equations 
admit a separation of variables. As a consequence, they admit simultaneous eigenfunctions. Consider the 
set of n Schrodinger equations for a state ^(a\f3): 

<j k {B)- T(P'\A)^(a\p)=T k {P'\A)^(a\p), for k = l,---,n (56) 

with eigenvalues <Jk(B). (Recall that P' = (— l) n P). As shown in Appendix E, we can look for eigenfunctions 
^(al/?) in a factorized form: 

*(a|/3) = ^(a 1 |/?Ma 2 |/3)---^K|/3) (57) 

provided the function ?/;(a|/3), which depends only on one of the a variables, is solution of the following 
separated difference equation: 

™ / B — A \ 

P 3 ■ ${ aj |/3) = + iC\l3) = II { Bk k +qAj l/3) (58) 
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This is the quantum analogue of the classical separation of variables. A remarkable fact of eq. (pq) is that 
its solution can again be factorized into product of functions depending separately on only one Pk- 

n 

fe|j8) = 11^1/3*) 
fe=i 

where ip(a\/3) satisfies eq.(|33|). This probably reflects the duality symmetry between the a and (3 variables. It 
is clear that the function ip(a\[3) is defined by the difference equation up to multiplication by any if -periodic 
function. The way to fix this ambiguity is the following: one requires that the function ip(a\/3) is regular for 
< Im(a) < 2ir, ip(a\{3) = 0(exp((7r/f — l)a) when a — > +oo. 



5.4 Exact form factors in the reflect ionless case. 

Let us return to the case £ = tt/u for integer v. In this case the wave-functions ^(a) are 27ri-periodic, that 
is why we denote them by *&(A). Moreover ^(A) are polynomials of A, so we take the Hilbert space as the 
space of polynomials. The operators Pj are identically equal to 1 because they correspond to shift of a's by 
2ni. Thus the formula for the scalar product must be of the form 

(*i|*a)= / dA x --- f dA n ]](A 2 ~A 2 ) 9\(A)Q( ai , ■ ■ ■ a„)* 2 (A) 
Jc Jc iK j 

In the reflectionless case the contour C are small contour around the origin. We recall that for f = ^ with 
v integer the function ijj{a\f3) is: 

v-l 

i>(a\P) = l[(B-Aq-i) 
i=i 

It is time to discuss the local operators in this A-representation. As we have seen from the exact quantum 
formulae, in order to insert the primary field $ TO into the matrix element one has to put under the integral 
the expression 

i\^f{A T f{B^ 

On the other hand, in the classical theory the fields $ m = e %mv are represented on the n-soliton solutions 
by = Ilj Ilj BJ™- After the quantization we have a self-adjoint operator H n = H n (A, P) such that 
H n ^(A, B) — (Y[j Bj)^(A. B). So, the comparison of the quantum and classical formulae shows that the 
classical expression Yl A r j L H~ m must be ordered for quantization as follows 

$ m (A, P) = H- m/2 {A, P) |j] Af \ H~ m ' 2 (A, P) (59) 

This ordering prescription ensures that & m is a real field, = $ m , since H n and Aj are hermitic. We 
hope that the same notation for the operator $> m acting in all the space of states of RSG or restricted to 
n-soliton subspace is not misleading. 

Now we can fix the function Q(a). The contours of integration are drawn around A = 0, so Q(a) has to 
be taken as Y[ a k ™ fc ■ positive powers of ak would give zero scalar product. It is also clear that the m^s must 
all be different for the anti-symmetry coming from J^[ i<; • (A? — A 2 ). Considering the form of the function tp 
one realizes that if one of the rrik is greater than n then the matrix element corresponding to the operator 
$i vanishes because the contour of integration with respect to Ak can be moved to infinity. Thus we are left 
with the only possible choice 

m k = k 

This is exactly what we have in the formulae known from bootstrap. 
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6 Concluding remarks 



Let us describe possible directions of future developments. 

We have constructed the local integrals of motion in terms of the operators A and P. One can consider 
the descendants of the primary fields with respect to these operators. However, to consider the full space of 
local fields we need to construct the Virasoro algebras in terms of A and P and to consider the descendants 
of the primary fields. We are sure that it can be done. In this way we must be able to identify the SG local 
operators described in p3] with those coming from the CFT description. For these computations we do not 
really need to study deeper the situation of generic coupling constant since the formulae from the subsection 
5.2 are absolutely general. 

There is another point for which the consideration of generic coupling constant is important. We have 
seen that there is a dual Weyl algebra composed of a and p which commute with the operators A and P. It 
is easy to argue that the non-local integrals can be expressed in terms of a and p. The commutativity of local 
and non-local integrals follows from the commutativity of A, P with a, p. The non-local charges represent 
the quantum loop algebra U-^sl^)- It would be very interesting to find their expressions in terms of a and p. 
On the other hand the non-local charges in the conformal limit correspond to screening operators. 

We want also to remind the reader that we were actually not working with complete form factors. We 
omitted certain multipliers which are the same for all operators, the normalization of the wave functions in 
the logic of this paper. We hope to explain how this piece appears in our approach in further publication. 
Here we just would like to stress the analogy with the method of orbits in coadjoint representations. Every 
orbit is quantized independently giving an irreducible representation of the group, but combining them in 
the regular representation requires to take into account the Plancherel measure: this is exactly the analog 
of the omitted normalization factors. 

There are amusing coincidences between many tools used in this paper and those existing in the works 
on lattice quantization of SG |2jj] : the Weyl algebras, the functions of the type of ip(a,(3). The latter are 
now called quantum dilogarithms because they provide a deformation of dilogarithm functions. It might be 
possible that this coincidence is not occasional. The classical soliton is similar to the step-function, so one 
can imagine that exact quantum mechanics of n-solitons is related to the theory on the lattice with n sites. 

These is a also a no-empty intersection with recent works on the Calogero models, Mac Donald poly- 
nomials and affinc Hccke algebras. We realize that techniques very similar to those we used to deal with 
the quantum r-function, i.e. the generating function of Mac Donald difference operators, may be used to 
separate the variables in these operators. One potential application of this remark would be a more detailed 
description of the algebra of the Mac Donald polynomials. 
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7 Appendix A: The analytical variables and finite zone solutions. 



Let us explain the origin of the A variables as a remnant of the parametrization of the finite zone solutions, 
when they degenerate to the soliton case. We shall do it in the simplest case of the KdV equation, cf eg. 
p7| . We start with an hyperelliptic curve T of genus n, and a divisor D of degree n on it. 

2n 

T:s 2 = R(X), R(\) = H(\-\ j ) 
D : = (ui,u 2 ,- ■ ■ ,v n ) 

We describe T as a two sheeted cover of the A plane. We put cuts on the real axis on the intervals (—00, Ao) 
and (Aaj-i, X21), i = 1, - ■ ■ ,n. The quantities z/; in the divisor D denote the A coordinates of the points of 
the divisor. One should keep in mind that to specify the points themselves, one has to choose the sheet 
above X = v%. With these data we construct the Baker- Akhiezer function which is the unique function with 
the following analytical properties: 

• It has an essential singularity at the point P+ above infinity: ^(^i A) = e kx (l + 0(1/ k)) with k = \fX. 

• It has n simple poles outside P + . The divisor of these poles is D. 

Considering the quantity ~d x ip + A^, we see that it has the same analytical properties as ip itself, apart for 
the first normalization condition. Hence, because ip is unique, there exists a function u(x) such that 

- d 2 x ip + u(x)ip + A^> = (60) 

We recognize the usual linear system associated to the KdV equation. One can give various explicit con- 
structions of the Baker- Akhiezer function. The most popular one is in terms of theta functions. However, 
for our purpose, another representation is more suitable. Let us introduce the divisor Z(x) of the zeroes of 
the Baker-Akhiezer function. It is of degree n: 

Z(x) := Oi(ar),M 2 (», • ■ ■ ,MnO)) 

One can find the equations of motion for the divisor Z(x). Consider the function d x ip/tp. It is a meromorphic 
function on T, it has poles at the points Hi{x) and behaves like k + 0(l/k) in the vicinity of the point P+. 
Hence we can write 



d x jj _ y/RjX) + Q(x,X) 

* nLi(A-^w) [ ' 

where Q is a polynomial of degree n — 1 in A. We determine Q by requiring that has a pole above 
A = fii(x) only on one of the two sheets (say y/R(X)). Then 



Q(x,ih(x))=t/R(jjh(x)) (62) 

Hence 

Q(x, A) = ^g i^ - (63 ) 



On the other side, in the vicinity of fJ>i(x), we have: 

d x ip d x fi t (x) 



ij) X- (m(x) 
Comparing eq.(J6l|) and eq.(pj), we get the equation of motions: 



0(1) (64) 



0^( X ) = -2 n f?f» ( „ (65) 
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One can now reconstruct the Baker- Akhiezer function itself. Indeed, inserting eq.(|6^) into ( |6l| ) we get: 

1 



ip P(X,x) 2^X-fii(x) 



where the polynomial P(X,x) is defined as P(X,x) = Yii(^ ~ Therefore |27 

i/j(X,x) 



ip{X,x ) 



1 P(\x) 
P(X,x ) 



exp 



dx 



(66) 



One can also reconstruct the potential u(x) directly in terms of the data /ii(x) and A^. Inserting back 
eq(|66|) into eq. (|60|) , we get the polynomial identity 

R = -]-PP" + -P' 2 + (u + X)P 2 
2 4 



Comparing the terms A 2 ™ we obtain: 



2» 



i=l i=0 

We are now ready to analyze the soliton limit. It corresponds to the following limiting configuration: 



A = 0, A 



2.7-1 



A 



2j 



m 



X = A 2 ; = A? 



(67) 



(68) 



In this limiting configuration, eq.(|67|) reduces to eq. (pl|) . Under the condition (|68|), the full curve T becomes 
the A plane but with the points Bj and — Bj, j = 1, • • • , n identified. In this limit we have: 



\=A 2 



= -A\{{A 2 -B 2 ), P(X,x) 



\=A 2 



A 2 



3 = 1 



The equations of motion for the Ai become (compare with eq.(20)) 

d^ = l[(A2-B*)l[^-p 

We can also obtain the degeneration of the Baker- Akhiezer function. Noticing that 
y/RW . 1 1 



(69) 



P(X,x) 

we can integrate eq.(|6^) to get 

ip(A, x) 



-A 



X=A 2 



9 51 ( A - 



A- A, A + At 



9 x Ai 



n 



ip(A,x ) f = \\A- Aj(x Q ) J 



A- A 



,(xy\ 



exp (— A(x — xo)) 



(70) 



As we see, in this limit the Baker- Akhiezer function contains a single exponential factor instead of two as 
one would expect from the second order linear equation eq.(|60[). This correspond exactly to the fact that 
for soliton solutions the potential u in eq.(p0|) is reflectionless. It is manifest in eq.(70) that the variables Aj 
are the zeroes of the Baker- Akhiezer function. 

As we said, the points Bj and — Bj are identified. This means that the Baker- Akhiezer function should 
satisfy the conditions ip(Bi, x) = ?/>(— i?j, x), i = 1, ■ • • n. Writing theses conditions, we get 



— -Z jo 



n Bi + Ajix) 



(71) 



where Yj depends only on the initial conditions at cc = xq. Comparing with eqs.(|]J^), we recognize the 
solution of the equations of motion. 
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In the sine-Gordon soliton case, we identify x = x_. The x + dependence is simply reintroduced by 
replacing in eq.(f7l|) 

e 2BiX _^ g 2B i i_+2B: 1 i + 

Let us discuss now the motion of the divisor Z(x-,x+). In the finite zone case, the point fj,i(x-) has a 
quasi-periodic motion on the real axis in the interval \2i-2 < fJ-i(x-) < ^2i-i- 

Consider the x+ motion in the soliton case when we have only two points. When x + = —00, the points 
A\ and A2 start from B\ and B2 respectively. When x + increases, the points Ai and A2 start to move to 
the right. Notice that the point A — 00 is regular for eq.(|69|), so the point A2 passes smoothly from +00 
to —00, and then continues to move towards — B2 which is reached at some time x*. At the same time x#, 
A\ reaches B2 so that in the right hand side of eq.(|7l|) the pole at A2 = — £?2 is cancelled by the zero at 
A\ = B%. Hence everything remains finite for a finite time x*. At this time, the point A\ jumps to — B2 and 
ends its motion at the point —B\. Similarly, A2 jumps to B2 and continues its motion up to — B2 again. 
The case of generic n is similar. Altogether, Aj starts at J3j and ends at —B. L . 

Eq. (ffoj) is easily related to the Jost solution of eq.(]60|). In our context, the Jost solution is defined by 
the normalization condition lim^ _^_ oc ^)j ost {x-) = cxp(— Ax J). Hence, we find 

x ^-oo IpBaker(XQ) f ^ \ A - Bj J 

where we used the fact that lim Xo ^- 00 Aj(xo) = Bj. Therefore Bj are the poles of the Jost solution and 
Aj(x) are its zeroes. 

Finally, we would like to discuss the important question of the Poisson structures of the KdV equation. As 
we know, we have a whole hierarchy of these structures. One can describe the restrictions of the symplectic 
forms to the manifold of finite zone solutions p^J. Let fi( fe ) be the restricted A:-th symplectic form. Then we 
have in terms of analytical variables p9j, ^(J : 

nW=£dP0*)A-T& W 

i=l Pi 
where V{\) is the pseudo momentum defined by: 

ip{X,x = L) 



V(X) = log 



ip(\,x = -L) 



To compute the quasi-momentum in the soliton limit, we choose the normalization point to be xq = —L, 
and we send L — > 00. According to the previous discussion, we have Ai(L) — > — Ai{— L) — * Bj. Using 
eq-©, we get 



V(A) = -2LA - log J] mod 11 

3 

Hence, 

n 

ow= 2 ^di og n 



i=l 



A 



Bj+Ai) Aj k ~ 3 



The form used in eq.(p5|) corresponds to k — 2 i.e., to the second Hamiltonian structure of KdV. We recall 
once more that the second Hamiltonian structure of the KdV equation is precisely the Virasoro algebra. 



8 Appendix B: From the {X, B} to the {A, B} variables. 

Before explaining the proof of the formula for the r-functions in the {A, B} variables we need to gather a 
few facts concerning the r-functions. The r-functions are defined by the determinant (|^): r± = det(l ± V). 
Most of the proof will be recursive using the recursion relation (||) satisfied by them, cf eg. ref. j2(| : 

r^(X\B) = T^iXlB) +T^{pl n {B)X k \B) X n . (73) 



22 



For comparison with the quantum formula, it is also useful to know the explicit expression of the r-function 
not in the Xj variables but in the variables Yj — Xj Ylk& f b^b^ ) '■ 

n 

t^(y\b)=i+j2 e iw^-n* (74) 



o=l rc{i, •••,«} iei ieJ 
/|=p jg* 



with " " 



U — Bi + Bj ■ 

Let us now prove the formula ( |l2] ) for the r-functions. We recall them to ease the reading: 

rj n) = 2" I I I Bj I X "" J ' - , J/± "^.\ — - ( To ) 



r («) = 




The upper index n refers to the n-soliton solutions. Eqs.(75) are two identities between rational functions 



of A and B once the expression of X as a function of A and B 

An) _ v (n) -Q / Bj - B k \ _ JT ^ -Sj - ^ 

fe#j + ^ fc / fc=l 



has been inserted into the r-functions. 

Let us denote by t^\a\B) the r-functions with Xj expressed in terms of A and B. Since t^\x\B) are 

symmetric in Xj, so is t± (A\B) as a function of Aj. Since permuting the B permutes the X, the functions 

t±^ (A\B) are also symmetric in the B. Thus, the identities ( |75| ) are equalities between rational functions 
symmetric in A and B. 

Let us first show that t± (A\B) has poles only at Aj + B k = 0. In view of the explicit expressions of the 

r-functions and of the Xj , t±\a\B) has potentially simple poles at Aj + B k = and Bj ± B k = 0. The 
expression (74) of the r-function in terms of Y shows that there are no poles at Bj + B k = 0. Similarly, the 
expression (W of the r-functions in terms of X shows that the potential poles at Bj — B k = are associated 
to Xj and X k . Using twice the recursion relation ([73]) shows that these poles cancel against 0fj(B). Thus, 

t± (A\E) can be written as 

Q ( ±\a\b) 



>{A\B) 



U hk (Aj+B k ) 



where Q±\A\B) are polynomials, symmetric in A and B, and of degree at most n in each variables. 

To prove the identities ( JT5| ) it is then enough by symmetry to show that the functions have identical 
residues and that they coincides at particular points. Again by symmetry, it is enough to check it for the 
pole at A n + B n = and at the point A n = B n . To keep the size of this Appendix reasonable we shall 
describe the proof for r + only. The proof for r_ is similar. We shall prove it by induction assuming that the 
identities are true up to n — 1 solitons (they are obviously true in the one soliton case). The recursive proof 
is based on the three following relations: 

i) The r-function satisfies the recursion relation ([73]) . 

ii) The r.h.s of eq.(|7^) satisfies the following recursion relation: 

rhs (n) = ( 2B n \ TT (An + Aj)(B n + Bj) ( „_ x) 

{A n + Bjj.J k (A n + Bj)(B n + A j ) (77j 

iii) As functions of A and B the variables Xj satisfy the following recursion relations: 

,(„) _ / B n + Bj \ ( Bj - A n \ (n-l) 



X ' - ){tTTj X > /«■ i - V ■ ■ ,» - 1 (78) 
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where xj"" 1 -* is independent of A n ,B n . 

Let us first compare the values of both side of eq.((75j) at the point A n = B n . Using eq.(f77|), we have for 
the r.h.s.: 



rhs {n) 



A n =B„ 



For the other side we remark that Xn vanishes for A n = B n . Therefore the recursion relation ( |73| ) for the 
t- functions implies that: 



lhs {n) 



lhs {n - 



A n =B„ 



Thus, both sides of eq. (|7q) coincide at A n — B n . 

Let us now compute the residue at the simple pole A r , 
pole of Xn a t A n = —B n . Its residue is: 



B n . The pole in t^\A\B) comes from the 



Res(X^) 



A n =-B„ 



(B n 




(B n 






-bJ 


\B n 





Furthermore, eq.([7q) gives for 1 < j < n— 1: 



X 



(n) 



A„=-B, 



= P7n(B)X 



(n-l) 



Therefore the factor /3? (B) cancels in the recursion relation ([73|) and we get: 



Res(T^\A\B)) 



On the other hand, eq.(|77|) implies: 

Res{rhs {n) ) 



A n =-B n 



Res(X^) 



>(n-l) 



A n =-B„ 



(A\B)) 



A n =-B n 



2B n Y[ 

Res(X^) 



(A 3 - B n ){B n + Bj) 



(Bj - B n )(B n + Aj) 



■ rhs (n - 1) 



\A n =-B n 



Thus, the residues of both sides of eq.(75) at A n = —B n coincide by the induction hypothesis. This concludes 
the proof of the identities ([75|). 

Let us now consider the formula expressing the symmetric functions <Jk{A) as functions of the {X, B} 
variables. As pointed out in the main text, the defining relations (|^) or ( f76| ) can be considered as a system 
of equations for the symmetric functions of A. We claim that the solution to this system can be written as: 



(Tk{A) 



T k (X'\B) 
r(X'\B) 



(80) 



where X' = {-l) n X and 



Tk(Xi, ■ ■ ■ , X n \Bi, ■ ■ ■ , B n ) — 

Bi x B l2 ■ ■ ■ B lk r(Xi, ■ ■ ■ , — Xi 1} ■ ■ ■ ,—X tkl ■ ■ ■ ,X n \Bi, 



. B n ) 



il<i2<---<ik 



Let us compare eqs.(^) and (|l4|). A quick look reveals that they are identical provided we exchange A with 
B and Y with P. Hence, the proof of eq.(|80|) is identical to the proof of eq. (|l5|) , which is a limiting case of 
the quantum formula whose proof is given in Appendix E. 

This somewhat mysterious formula can be related to well known results connecting the Baker function and 
the tau function. One first reintroduce all the times by the substitution X{ — * Xi(t) = Xi(0) exp(2^(Bj, t)), 
where £(A, t) = J2i A 2l ~ 1 t 2 i-i- Next, using the generating function eq.(p7|) we have 



uB-j 



r(X\B) 
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On the other hand, the multi time Baker function reads 



VKA t) = yrf 1-Mt)/A \ e(A , f)+C(i i,tW) 

i>(A,t(°)) *-l\l-Ai(j0))/A) 



Combining these two formulae, we can identify the Baker function as 



t(X\B) 

Since 1^5773 = ex P ( 2 E* 2^=1 (t) 2 " -1 ) we nnd 

= T{t2n - 1 + T ^^ =r i e -«^) 
and we recognize the well known Sato formula. 

9 Appendix C: Information about the integral formulae. 

In this appendix we explain how the formulae for the form factors given in the Section 3 agrees with the 
conventional ones We have to show how to reduce the number of integrals and the degree of polynomials. 

In order to reduce the number of integrals by one let us consider the integral with respect to A n in which 
we would like to move the integration contour to infinity. When A n — > 00 one has 



So, the integral over A n is 



g ^(A^Bj) ]J(4 2 - A 2 n )a- n = A- 2 ( 1 - i±|A- 1 £ B 3 + 0{A; 2 ) ) (81) 



, A n -!,A n \B B 2n ) dA n (82) 



It is important that there are no contributions with a n = A 2 ^ which means that the integral is essentially 
independent of the coupling constant (it depends on v only through the constants like ^-j- in (|8^)). The 
first two terms written in ( |S2"| ) are sufficient to calculate this integral for 

L 9l {A u ■ ■ ■ , A n \B u • • • , Ban) = II A i II B 3 ' 

3 j 

and 

{A lr --,A n \B u ---, B 2n ) = Y[ A) H Bj 1 

3 3 

the results being respectively 

n—1 n— 1 

«=i 3 3 »=i j 

This calculation gives agreement with the formulae from ||. To calculate the integral over A n for higher 
operator $ m one has to take into account higher contributions to ( |Sl| ) , it would be nice to find these primary 
operators among those described in . 

Let us now explain why the reduction of the degree of polynomials is possible. The formula for /© is 
composed of one-fold integrals of the type 

2n 

Y[^{A,B j ) a - k L(A)dA 

a i =1 
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One can reduce the degree of the polynomial L(A) using the following circumstance. Due to the fact that 
the function ip{A, B) satisfies the equation ( |33| ) 

^(Aq,B) = ^-^{A,B) 

the polynomials of the form 

L(A) = M(A) \\{B 3 +A)- qM(Aq) \\{B 3 - A) ~ (83) 

3 3 

for any polynomial M (A). Here L(A) ~ means that L(A) produces zero when substituted into the integral. 
This allows for any given polynomial L(A) to find such a polynomial L'(A) such that L(A) ~ L'(A) and 
deg(L') <2n—l. To do that one has to find a polynomial M(A) such that 

L(A) = M(A) HiBj +A)- qM(Aq) ]J(B, - A) + L'(A) 

3 3 

with deg(L') < 2n — 1 . This is always possible by induction. 

10 Appendix D: A semi-classical limit. 

We start with 

v—l 

i>(A,B) = l[(B-Aq-iy, q = e i % 

3 = 1 

We recall the obvious formula 

n— 1 n=0 

= ~f f(x)dx-~(f(0) + f(N)) + O(A) 

Using this formula we find 

logV = - / log(B-Ae- lx )dx--(log(B-A)+\og(B + A)) + 0(A) 
* Jo 2 

The integral defines an analytical function of A in the plane with a cut which is a semi circle from B to — B 
in the upper half plane. Remark that for A = 0, this integral is real and its value is 7rlogi?. To proceed, we 
perform some formal manipulations on the integral. We have 

f* C A dA P A dA f~ A dA 

/ log(S - Ae- ix )dx = i log(B + A)—=i \ og (B + A)—-i \og{B + A) — 
Jo J-A A J Q A J A 

A / B + A\ dA 



o 



where the last term has been added to normalize the function by its value at A — 0, thereby fixing the 
ambiguities of the formal manipulations. Putting everything together, we get 



B v 

ib(A, B) = ; exp 
VV ' 1 VB 2 - A 2 1 



v f A , {B-A\ dA 



l «Jo l0S \B + Aj A 



which is the result quoted in the text. 
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11 Appendix E: A proof of the quantum separation of variables. 

Before proving the separation of variables in the quantum theory, we present a proof for the formula of the 
generating function T'"'(u) = J2 k ^T^. ^ s ^ ^ s defined in (|52"|), the quantum r- function satisfies the 
following recursion relation: 

T™(P\A) =: T^( lkn P k \A k ) + (J] lnk )T^ (%^P k \A k ) P n : (84) 

where 

qAi + Aj _ q- x Ai+Aj 

7y ' = 4- 4- ' and ^ = A- A- = _97ii 

The double dots : : mean writing the P's on the right. In the classical limit q — > 1 this is equivalent to the 
relation (Sj. 

Let Tj* be the operators defined in eq.([5l|) in the n solitons case. By convention we set: 7jj ' = 
T^ n \P\A) and T fc (n) = for k < or k > n. From eq.(p4h, we deduce a recursion relation for the T k : 

, ~(n-l) / , . ~(n— 1) \ 

T W=f fe (n - 1) +f fe Pn + A.^ 1 '-^*-! PnJ (85) 

where 

_ _ ~(ra-l) ~(n-l) 

bumming up eqs. dH) and defining T^" 1 )^) = £ fc u k T k {n ~ 1] and T (u) = £ fc u k T k , we get: 
T (B) («) = (l + tMnjr^CuJ + Cl-guA.JT (u)P„ 



(?<-»(«) + ?- W P„(i^))( 1+ ^„, 



Comparing with the recursion relation ( p4[ ) satisfied by the quantum r-function 7~(P|A) proves the result 
quoted in eq.(|53|). 

Furthermore, the hermiticity properties of the quantum r-function t(P\A) implies: 

n ^~ ^ 1 — A 

T(u)t = ;Q(l + uA,)T(gP(u)|A), u»t/i P 3 -(«) = —- -/P; 

3 =l 

We can use the generating function T(u) do write the commutativity relation ( ^5|) in an alternative form: 

n ^ n ^ 

T(P(u)\A) ■ + uA 3 )(l + vAj) ■ T(qP(v)\A) = T(P(v)\A) ■ J]^ 1 + vA 3 )(l + uA 3 ) ■ T(qP(u)\A) (86) 

3=1 3=1 

Let us now give a proof of the quantum separation of variables. Assume that we are acting on the 
quantum hamiltonian H k with a wave function ty(a\/3), eq.(|57|), satisfying the difference equation (p8[). 
Since the momenta operators P, have been ordered to the right in the definition of the quantum r-functions, 
this wave function will be an eigenfunction with eigenvalue <J k {B) if the following relation is true: 

a k {B) r^(pW|A) =T k (n) {P {n ^\A) (87) 

with 

p(«) = / Aj - B k 



k=l 



qA 3 + B k 
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Recall that we define the Hamiltonians using P' = (— l) n P. In eq. (|S7|) we used the same notation for 
the quantum operator T^ n \P'\A) and the c-number function obtained by inserting the values ( |S8| ) of the 
momenta operators PL We recall that the functions T^ n '(P\A) are defined by: 



T k {n) {P\A) = Yl A il A i2 ---A ih r(P 1 ,---,-qP il ,---,-qP il 



■ P n \A) 



(89) 



%X<%-2 — <lh 



Once the specific values 



of the momenta have been plugged into eq.(|87|), this is an identity between 
two rational functions in the variables A and B. Both the l.h.s. and the r.h.s. of ( |S7| ) are symmetric functions 
in A and B. They only have simple poles at qAj + Bk = and Aj — Ak, the former comes from P^ and the 
latter from jj). . Due to their behavior at infinity, to prove eq. (p7[) it is enough to check that these rational 
functions have the same residues at their poles and that they are equal at specific points. By symmetry it 
is enough to check the residues at qA n + B n = and A n = A n -i, and the values at the point A n = B n . We 
shall do it by induction assuming that the identity ( |87| ) is true up to n — 1 (it is of course true for n = 1). 

i) Consider first the residue at qA n + B n = 0. This pole is associated to P^\ To compute the residues 
of both side of eq.([S7|) we may use the recursion relation (|S4|). Using the fact that: 



1 

J Jn j qA n +B n =0 

the recursion relation ( p7| ) implies: 
Res(lhs) 



P 



(n-l) 



for 



J 



1. 



, n — 1 



qA n +B n =0 



(<t<»-'>(B) + B„4':- 1 1, (B)) 

- ((-»)-' n(^)-^-"^"-' 



\A) Res(P^) 



Similarly applying the recursion relation to the rhs, but distinguishing whether A n is a marked point in the 
sum (R9) or not, we get: 



Res(rhs) 



qA„+B„=0 



Bn 



\qAj + Br, 



T (n-l) (pin -l )lA) 



Res(P^) 



Comparing these formula gives the equality of the residues by the induction hypothesis. 

ii) Consider the residues at A n — A n -\. They are both vanishing, and therefore equal. The proof of the 
vanishing of the rhs residue is similar to the proof of vanishing of the lhs residue, so we shall only present the 
latter. The potential pole at A n = A n -i is associated to j n ,n-i- Therefore, its residue may be computed by 
using twice the recursion relation d§4|). We get: 



Res(lhs) 



A n =A n - 

This vanishes since we have 



= const, ( T ^\J^P 3 \A)P n ^ T^('^zlp j \A)P n ) 

1 \ / i,n— 1 ^ in / 



p 



(«) 



= p 



(n) 



iii) Consider now the value of the functions at the point A n = B„ . At this point we have P^ = and 



P 



(n-l) 



A n =B„ 



for 



■ ,71 - 



Therefore, the recursion relation (p4[) yields to 



Iht 



A„=B„ 



= (a^- l \B)+B n a k n _- 1 1 >(B)) T^- X ){P^\A) 
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for the left hand side and, 
rhs 



= Z 



(n- 



for the right hand side. It clearly appears that they are equal by the induction hypothesis. 
Collecting the points i) to iii) proves the quantum separation of variables. 
The classical formula corresponds to the case q = 1. 
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